First-principles study of adhesion at Cu/Si02 interfaces 
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The structural, electronic, and adhesive properties of Cu/Si02 interfaces are investigated using first- 
principles density-functional theory within the local density approximation. Interfaces between fee 
Cu and a-cristobalite(OOl) surfaces with different surface stoichiometries are considered. Interfacial 
properties are found to be sensitive to the choice of the termination, and the oxygen density at the 
substrate surface is the most important factor infiuencing the strength of adhesion. For oxygen- 
rich interfaces, the O atoms at the interface substantially rearrange after the deposition of Cu 
layers, suggesting the formation of Cu-O bonds. Significant hybridization between Cu— d and O— p 
states is evident in site-projected density of states at the interface. As oxygen is systematically 
removed from the interface, less rearrangement is observed, implying weaker adhesion. Computed 
adhesion energies for each of the interfaces are found to refiect these observed structural and bonding 
trends, leading to the largest adhesion energy in the oxygen rich cases. The adhesion energy is also 
calculated between Cu and Si02 substrates terminated with hydroxyl groups, and adhesion of Cu to 
these substrates is found to be considerably reduced. This work supports the notion that Cu films 
can adhere well to hydroxyl-free Si02 substrates should oxygen be present in sufficient amounts at 
the interface. 



PACS: 68.35.-p, 73.20.-r, 68.47.Gh 
I. INTRODUCTION 

Adhesion of thin metal films to glass substrates have 
continued to be the subject of intense study for many 
years because of its importance to the large-area elec- 
tronics industry. Copper in particular, a noble metal 
with a high bulk thermal and electrical conductivities, 
and apparently low electromigration rate,^ is an excel- 
lent candidate for interconnects in integrated circuitry on 
dielectric substrates, and thus the formation of a strong 
and highly reliable interface between copper and glass has 
been an issue of primary importance. Unfortunately, Cu 
films are often reported to bind rather poorly to oxide- 
based glass substrates, and films that do adhere at 
room temperature can show diminished adhesivity after 
thermal cycling, an effect often attributed to differ- 
ences of thermal expansion between Cu and the oxide 
substrate. * In many studies, poor bonding is frequently 
explained by noting that since Cu does not reduce Si02, 
a graded oxide layer facilitating the adhesion is unable 
to form. Chemically, copper possesses a half-empty s 
shell and a reasonably tightly-bound and filled d-shell, 
and thus is less reactive than aluminum and many tran- 
sition metal films, ^ which oxidize well and are observed 
to adhere more reliably. Thus it is not surprising that 
adhesion is often improved experimentally by applying 
transition metal intermediary layers^'^ prior to Cu film 
deposition, or alloying the Cu films with small amounts 
of Mg or A1.4.5'9"i2 

Adhesion of copper films to oxide glasses is, of course, 
expected to be highly sensitive to the temperature, oxy- 
gen partial pressure, and the condition of the glass sur- 



face prior to deposition. A common feature of many stud- 
ies of Cu films on glass is that they were performed under 
conditions amenable to surface passivation either before 
or during deposition of the copper film. For example, Cu 
films sputtered under high purity conditions in a vacuum 
do adhere to glass substrates. ^'^'^ In particular, Ohmi et 
al7 demonstrated robust adhesion of copper to glass af- 
ter removing hydroxyl (OH) groups from the interface in 
situ just prior to deposition. This suggests that dangling 
bonds at moisture-free glass surfaces would be readily 
saturated by Cu-O bonds. 

Glasses important to industry can be chemically com- 
plex, often containing a substantial amount (perhaps 
~ 30%) of aluminum, boron, and other alkaline earth 
elements. Their primary component, however, is al- 
most always Si02, which by itself forms one of the sim- 
plest glasses, amorphous silica (a-Si02). In the interest 
of simplifying the chemistry, we shall restrict our focus 
to chemically-pure Si02 glasses for the duration of this 
study. Amorphous silica has no long range order yet re- 
markably, on length scales comparable to a Si-0 bond, 
it is nearly perfectly ordered: the bond-length does not 
vary appreciably from 1.61 A; each Si atom is tetrahe- 
drally coordinated with O atoms, and the bonds are pri- 
marily of a covalent nature. At length scales between 
5 A and 8 A, minimal variations in 0-Si-O and Si-O-Si 
bond angles result in an appreciable degree of interme- 
diate range orientational disorder. Previous work using 
molecular dynamics has been performed on both bulk vit- 
reous silica and silica surfaces. ^^"^^ Much of it was based 
on semi-empirical potentials, which though able to 
simulate systems large enough to correctly capture the 



1 



effects of the intermediate range orientational disorder, 
are expected to be far less accurate when local condi- 
tions deviate from those in the bulk, such as they do at 
surfaces or interfaces. 

Since microscopically adhesion is related to the 
strength of the electronic bond between atoms at the in- 
terface, the local electronic structure at the interface will 
play an important role in understanding the reactivity 
of the metal with the oxide substrate. A first princi- 
ples treatment would therefore be most appropriate for 
examining the nature of local bonding at an ideal inter- 
face between the two systems, and indeed this approach 
has proved insightful in the past for studying adhesion 
between metal films and insulating substrates. Pre- 
vious first-principles work on the a— quartz surface^° ob- 
served significant changes in bonding near the surface; to 
our knowledge, however, the effects of metallic overlay- 
ers on Si02 surface geometries and electronic properties 
have yet to be examined from first principles. 

In this article, we assess the degree to which the chem- 
ical bonding at atomically-sharp interfaces reflects the 
empirical trends in the adhesion of metallic contacts to 
glass substrates through study of a simplified system, Cu 
monolayers on a— cristobalite, a crystalline Si02 poly- 
morph which although ordered does possess a density 
quite close to that of a-Si02.^^ Adhesion is studied as a 
function of oxygen surface coverage from first principles, 
using density functional theory within the local density 
approximation (LDA).^^ Our methods allow us to elu- 
cidate the role of local chemistry in binding the film to 
glass, but require us to neglect its inherent long-range dis- 
order and also, to a limited extent, its intermediate-range 
orientational disorder. Since we are chiefly interested in 
the local bonding properties, which are expected to be 
the same for both amorphous and crystalline silica, the 
loss of disorder is acceptable. We summarize the method 
and approximations used here in Sec. II; in Sec. Ill we 
describe the relaxed 1x1 Q;-cristobalite(001) surfaces con- 
sidered here, which include both stoichiometric and non- 
stoichiometric surface terminations, chosen so as to re- 
flect different possible oxygen coverages that may result 
from different deposition conditions. In Sec. IV, we de- 
tail results of atomic relaxation of Cu/Si02 interfaces, 
where the adhesion is found to be critically sensitive to 
oxygen density at the substrate surface, and where rela- 
tively strong bonding is reported for oxygen-rich surfaces. 
The effect of hydroxyl groups as a surface passivator is 
also briefly considered by calculating adhesion energies 
between copper monolayers and a hydroxyl-terminated 
Si02 surface. After discussing our results and their im- 
plications, we provide concluding remarks in Sec. V. 



II. METHODOLOGY 
A. Computational details 

Our flrst-principles density-functional calculations are 
carried out within the LDA, using the correlation en- 
ergy of Ceperley and Alder, as implemented within 
the Vienna ab initio Simulation Program^^""^" (VASP). 
We use a plane-wave basis set with a 29 Ry cut-off, 
and ultrasoft pscudopotentials'^^''^^ for Cu and H, projec- 
tor augmented-wave potentials'^^''^'* for O, and a norm- 
conserving pseudopotentiaP^ for Si. The Monkhorst- 
Pack fc— point mesh is chosen so that = 3x3x1 
for the 5Ax5Ax30A supercell which is used in this 
work. Also, we make use of a Fermi-distribution smear- 
ing with the temperature of ksT ~ 0.2 eV to facilitate 
rapid convergence. These computational conditions pro- 
vide good convergence of the structure, density of states, 
and energy differences given here, and are used through- 
out this work unless otherwise stated. 

It should be mentioned that the fe— point sampling dis- 
cussed above is not entirely sufficient for study of the 
bulk properties of copper. In fact, the total energy of 
bulk copper obtained with this sampling is converged to 
about 0.02 eV per Cu atom. However, a sparser mesh 
turns out to be adequate for this study since energy dif- 
ferences are converged. To demonstrate, we increased 
our sampling to rifc =4x4x1 for the supercell men- 
tioned above, recalculated total energies and forces, and 
found that the resulting difference in adhesion energy is 
less than 4% even in the largest case. 

B. Bulk properties of constituent systems 

Copper crystallizes in the fee structure under normal 
conditions and its lattice constant is 3.615 A at 291 
K.'^^ The LDA underestimates this lattice constant as 
is well known, and we obtain 3.532 A for the equilibrium 
value as shown in Table I, in agreement with a previous 
LDA study.^^ At equihbrium we find a = 4.975 A and c 
= 6.907 A for the tetragonal cell of bulk a-cristobalite 
(the _P4i2i2 structure), in excellent agreement with 
experiment"^* and previous LDA calculations. ^"'^^ Fur- 
ther, the computed Si-0 bond length (~ 1.606 A) is re- 
markably close to experiment 1.603 A), and the cal- 
culated 0-Si-O bond angle is near the tetrahedral ideal 
value, also in-line with measurement.^^ Conveniently, the 
choice of a-cristobalite circumvents spurious lattice mis- 
match when considering the interfaces with Cu overlay- 
ers: Compare the lattice constant a of a— cristobalite 
with y/2a of Cu in Table I. 
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TABLE I. Lattice constants of Cu and a-cristobalite ob- 
served in experiment and LDA calculation.''^ Wyckoff posi- 
tions, two Si-0 bond lengths, and four O-Si-O bond angles are 
also shown for a-cristobalite; the Si Wyckoff position (4a) is 
shown as (u, u, 0) and the O Wyckoff position (8b) as {x, y, z). 
The experimental data are taken from Refs. 36 and 38 for cop- 
per and a-cristobalite, respectively. 
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III. (001) SURFACES OF a-CRISTOBALITE SIO2 

As a starting point, we consider an a— cristobalite slab 
consisting of 5 Si02 layers, the thickness of which is 
^--^8.3 A. Translational symmetry is removed along the 
c-axis, defining two different (001) surfaces perpendicu- 
lar to this axis. We embed this a-cristobalite slab, along 
with about 20 A of vacuum, within a supercell having 
dimensions of 5 Ax 5 Ax 30 A. Oxygen number at the 
top of our Si02 slab is then changed to mimic different 
substrate terminations. Two nonbridging oxygen atoms 
terminate the bottom of the slab, and two additional H 
atoms are attached to the oxygens to remove dangling 
bonds. Forces on the atoms are computed using the 
Hellmann-Feynman theorem, and the positions of atoms 
arc then updated until the total energy reaches a min- 
imum. The oxygen atoms at the bottom surface, ter- 
minated by hydrogens, are fixed at their bulk positions 
in order to reduce size effects owing to the finite thick- 
ness of the Si02 slab, though the positions of each of the 
hydrogen atoms are permitted to relax. Previous first- 
principles calculations of a-quartz surfaces'^" and model 
Si/Si02 interfaces'*^ have indicated that about 5 A away 
from the interface, the local structural and electronic 
properties of the slab are bulk-like, and thus the size 
of our slab is sufficient to approximate both surface and 
bulk features. 

Figure 1 presents views of the relaxed 
a-cristobalite(OOl) surfaces. Supercells with slabs hav- 
ing three different surface terminations are considered in 
this work, which we label as Si-, 0-, and 00-terminated. 
The Si-terminated slab contains five Si02 units together 




Si-terminated O-terminated 00-terminated 



FIG. 1. Relaxed structures of the a-cristobalite(OOl) sur- 
faces. 

with two bottom hydrogens; the O-terminated slab is 
identical to the Si-terminated except that the dangling 
surface oxygen atom is added; the 00-terminated slab is 
again a modification of the Si-terminated slab in which 
two oxygens are added to the surface. Bond lengths and 
angles are not expected to remain bulk-like at the surface, 
where the choice of termination can leave non-bridging 
oxygens with dangling bonds at the surface (e.g., in the 
00-terminated cases). Indeed, we find that relaxation 
significantly changes the Si-0 bond lengths within two 
atomic layers from the surface for all terminations consid- 
ered, and the nature of these changes is strongly depen- 
dent on the termination. In the Si-terminated case, the 
additional charge carried by the dangling Si bonds forces 
the Si surface atom above the plane of the surface; as a 
result the Si-0 bonds supporting it from below elongate 
slightly to 1.64 A. In the O-terminated case, the surface 
O atom, bonded to only one other Si, relaxes downward 
into the surface, and its bond length decreases by more 
than 5% from 1.6 A to 1.5 A. In the 00-terminated case, 
in some contrast to the Si- and O-terminated slabs, the 
Si-0 bond lengths remain essentially the same as bulk. 
But the distance between O neighbors on the surface is 
also 1.6 A, which is very small compared with typical bulk 
O second neighbor distances of 2.7 A. Thus O atoms may 
evidently saturate dangling bonds through bonding with 
neighboring O atoms at this surface. 

As expected, the total cohesive energies of the super- 
cells decrease substantially with each additional O atom: 
adding one oxygen atom or two oxygen atoms lowers 
the cohesive energy by ~7.1 eV or ~11.3 eV per surface 
Si, respectively. Since the energy required to break the 
O2 bond is nearly 5 eV per molecule, this energy com- 
parison suggests that the 00-terminated surface may 
be the most stable of the three terminations, and the 
O-terminated surface the second. Evidently our first- 
principles calculations are consistent with the intuition 
that glass surfaces should be oxygen-rich. 
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FIG. 2. Relaxed structures of the Cu/a-cristobalite(001) 
interfaces. 



IV. CU/SIO2 INTERFACES: RESULTS AND 
DISCUSSION 

A. Interfacial structure 

We now add five Cu layers to each of the fuUy-relaxed 
(001) a-cristobahte surfaces, obtained as discussed above 
in Sec. III. The Cu layers are initially positioned on the 
surface so that a Cu atom in the lowest monolayer lies on 
the former symmetry axis of the a-cristobalite slab. The 
length of the supercell lattice vector normal to the surface 
(the c-axis) is fixed to 30 A as in Sec. Ill, a value we find 
to be large enough to accommodate up to five additional 
Cu layers while keeping the interactions between the sur- 
faces in neighboring supercells minimal. Each monolayer 
in our supercell contains four Cu atoms, and before re- 
laxation each Cu atom has four intralayer neighbors at 
a distance of 2.5 A. Fortunately, the lattice parameter of 
the 1x1 surface of a-cristobalite is well matched to that 
of the \/2 X \/2 surface of copper as shown in Table I (the 
mismatch is less than 1%). The small lattice mismatch 
is artifact of our approximate treatment of amorphous 
silica as crystalline. In reality, stress induced by lattice- 
matching is relieved through formation of defects and/or 
dislocations. (In the case of a truly amorphous substrate, 
stress may also be overcome through surface reconstruc- 
tion, a complexity we neglect in the present work.) 

The positions of all atoms in the supercell are again 
relaxed except for the bottom O atoms, and Fig. 2 shows 
the optimized structures of the interfaces for each ter- 
mination. The most significant reconstruction at the 
interface is observed in the OO-terminated case. The 
0-0 bond observed at the free a-cristobalite surface in 
Fig. 1 is broken by the deposition of copper, and the 
0-Si-O angle is changed from 59° to 104° with the Si-0 
bond lengths kept almost constant. This large rearrange- 
ment would indicate bond formation between interfacial 
O and Cu atoms, suggesting good chemical adhesion be- 
tween the oxide substrate and the Cu film. Each interfa- 
cial O atom has two Cu neighbors, and the Cu-0 bond 
lengths are computed to be ~ 1.9 A, strikingly similar to 
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FIG. 4. Change in electron density {Sn{z)) (see text) as 
a function of the depth z from the substrate surface. The 
positions of each atom are designated by open circles for Si, 
solid circles for O, open squares for Cu, and open triangles 
for H. 



Cu-0 distances found in cuprates,^'^ where each copper 
atom is coordinated by four oxygens and considered to 
be in a Cu^+ state. Similar Si-0 and Cu-0 bond lengths 
are observed at the interface between O-terminated sur- 
faces and Cu monolayers, although the Si-O-Cu angles 
are smaller than in the OO-terminated case by ~ 15°, as 
can be seen in Fig. 2. The magnitude of the reconstruc- 
tion at the 0-terminated/Cu interface is less than that 
at the 00-terminated/Cu interface. 

The interface between the Si-terminated substrate and 
Cu is completely different from either of the oxygen 
terminated interfaces. In this case, the interfacial Si 
atom possesses four Cu neighbors, each with bond length 
~ 2.4 A, and there is negligible atomic rearrangement, 
implying more metallic-like bonding. 



B. Electronic properties 

To investigate further the bonding properties of each 
interface discussed above, we compute the site-projected 
local density of states (LDOS) for atoms near the in- 
terfaces, and compare them with those obtained for 
bulk-like atoms deeper inside the slabs; the LDOS are 
shown together for each interface in Fig. 3. In the OO- 
terminated case, significant hybridization is observed be- 
tween Cu— d and 0—p states just below the Fermi level. 
A key feature is the slight bump in the 0—p LDOS, ly- 
ing about 1 eV below the Fermi energy, which becomes 



4 



Si-terminated 0-terminated 00-terminated 




FIG. 3. Local density of states (LDOS) projected by site and angular momentum. The solid lines show the LDOS of atoms 
at (or near) the interfaces, and the dashed lines show those of atoms in the middle of the slabs (the LDOS of an approximate 
bulk material). When computing the LDOS, a gaussian smearing method is used with a = 0.2 eV. The sphere sizes of each 
atom are taken as 1.11 A for Si, 0.73 A for O, and 1.38 A for Cu. 



progressively smaller with decreasing O content at the 
interface, and nearly vanishes in the Si-terminated case. 
Evidently strong hybridization between Cu— d and 0—p 
states is responsible for the significant reconstruction at 
the 00-interface. It is worth pointing out here that as 
the number of O atoms at the interface is decreased, the 
Si— LDOS around the Fermi energy increases, that is 
the interfacial bonding takes on a more metaUic char- 
acter. This behavior results from the hybridization of 
Si— sp states with the itinerant Cu— s states around the 
Fermi energy. In addition, we observe the bandgap com- 
puted within the LDA (an underestimate of the true 
bandgap) opens and gradually approaches its bulk LDA 
value for LDOS in the middle of the Si02 slab, confir- 
mation that our slab exhibits approximate bulk behavior 
away from the interface. 

To examine the degree to which changes in structural 
and bonding properties are confined to the interface, and 
to observe the mixing between the electronic states of 
the Si02 substrate and the Cu layers in more detail, we 
calculate the density difference from the superposition, 
namely, 

{5n{z)) = {nw{z)) - [(nsiOsC^)) + {ncu{z))] , 

where {niF(z)), (^^3102 and (ncu(z)) are the densi- 
ties of the Cu/Si02 interface, the Si02 substrate, and the 
Cu layers, respectively, averaged over the xy plane (par- 
allel to the interface). To obtain (nsi02(^)) (('^Cu(2:))), 
we simply remove the Cu layers (Si02 substrate) from 
the fully relaxed interface in the supercell, and then re- 
calculate the electronic structure self-consistently keep- 
ing the atoms fixed. The quantity {6n{z)) thus indicates 
the change in electron density resulting from chemical 



TABLE II. Adhesion energy (ideal work of separation) W 
of the Cu/a-cristobalite(001) interfaces. 



W (J/m^) 


Si-terminated 


1.406 


O-terminated 


1.555 


OO-terminated 


3.805 



bonding between the Cu layers and the Si02 substrate. 
(To simplify the analysis, we are neglecting additional 
changes stemming from the atomic rearrangements at 
the free surfaces of Si02 and Cu.) The results of calcula- 
tions of the {dn{z)) appear in Fig. 4. The {6n{z)) in the 
OO-terminated case most prominently deviates from zero 
around the interface, reflecting significant charge trans- 
fer between the two slabs. The charge transfer in the O- 
and Si-terminated cases is comparatively much smaller. 
In particular, in the Si-terminated case, the depletion of 
the density in the vicinity of the lowest Cu layer is not so 
drastic, suggesting that the localized Cu— d states have 
less influence on the bonding around the interface. 



C. Quantitative analysis of adhesion 

To assess the adhesive strength of the interfaces, we 
have calculated the ideal work of separation, or adhesion 
energy, per unit area W from 

W = {Es,o, + Ecu - Eif)/A, 

where -Esi02 > -E'cu , and Eif are the energies of the iso- 
lated Si02 substrate, Cu layers, and Cu/Si02 interface in 
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FIG. 5. Adhesion energy as a function of the position of 
Cu layers in the xy plane calculated for five different initial 
positions and for each of the three terminations. Each point is 
obtained by translating Cu layers [initially positioned on the 
substrate as described in Sec. IV A ( Aa; = 0, Ay — 0)] in the 
xy plane by (Ax, Aj/) and then fully relaxed. The abscissa 
(which quantifies the magnitude of the shift of Cu layers) is 
scaled by the supercell lattice parameter o=5 A. 

the supercell, respectively, and A is the area. Physically, 
the adhesion energy W is the work per unit area required 
to separate the interface into the Cu layers and the Si02 
substrate within a microcanonical process, and it can be 
considered a measure of the strength of the adhesion. 
For the purposes of comparison, all energies are calcu- 
lated using supercells of the same size (5 Ax 5 Ax 30 A) 
independent of whether it contains Cu layers, any of the 
Si02 substrates, or Cu/Si02 interfaces. Table II lists 
the adhesion energy for each case. The adhesion en- 
ergy in the 00-terminated case turns out to be much 
larger than those in the O- and Si-terminated cases, 
and comparable to values previously computed for other 
metal-dielectric interfaces, such as Co/TiC (Ref. 21) and 
Nb/Al203 (Ref. 22). We note that the computed ener- 
gies for the 00-terminated case, ^ 4 J/m^, are consistent 
with values obtained experimentally by Kriese et al.^ for 
100 nm thick films of Cu on Si02 using an indentation 
technique. In summary, the magnitudes of adhesion en- 
ergies computed here for each of the three different inter- 
faces reflect the trends witnessed above in their structural 
and bonding properties. 

1. Local atomic structure 

At this point, it is meaningful to investigate further 
a drawback of using crystalline a— cristobalite to model 
the amorphous substrate. When Cu layers are placed on 
an amorphous Si02 substrate, different local atomic con- 



(1/4,0) (1/4,1/8) (1/4,1/4) 

FIG. 6. Structures at the 00-terminated interfaces ob- 
tained as a function of initial positions of Cu layers with re- 
spect to the substrate. The figures within brackets denote the 
shift of the position of Cu layers within the xy plane (see the 
caption of Fig. 5) in unit of 5 A (the in-plane supercell lattice 
parameter) . 

figurations are possible at the interfaces because of the 
non-periodicity of glass. A principal simplification of our 
model is the imposition of translational symmetry: we 
are unable to directly assess the difference between the 
adhesive properties of crystalline substrates and those of 
amorphous substrates. We are, however, able to investi- 
gate the sensitivity of our calculated adhesion energies to 
changes in local atomic structure. As a first step, we com- 
pute the influence of the position of Cu layers, relative to 
the substrate, on the adhesion energy by translating the 
Cu layers in the xy plane. 

As shown in Fig 5, the adhesion energy turns out rather 
insensitive to the initial position of the Cu layers so long 
as the termination type remains unchanged. This result 
suggests that the adhesive properties realized in these 
a— cristobalite substrates may well be carried over, to 
some extent, to those of glass substrates. The insensitiv- 
ity of adhesion energy to the position of Cu layers does 
not imply that the detailed structure at the interface is 
less important for adhesion. In fact, as shown in Fig. 6 
for the 00-terminated case, atoms at the interfaces do 
move significantly relative to the unshifted case. Since 
the surface atoms have a greater freedom to choose their 
position than those in bulk, the adhesion energy remains 
rather constant for various initial positions of Cu layers. 

2. Local oxygen surface coverage 

Thus far we have investigated the Si-, 0-, and 00- 
terminations independently. Depending on the oxygen 
partial pressure, however, these terminations are gener- 
ally expected to coexist on the substrate surface, inde- 
pendent of whether the substrate is crystalline or amor- 
phous. To examine the effects of oxygen surface densi- 
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Oxygen density (nm' ) 



FIG. 7. Adhesion energy as a function of oxygen density at 
the substrate surface.** The densities 0, 2, 4, 6, and 8 nm^^ 
correspond to the Si-, (Si,0)-, 0-, (0,00)-, OO-terminated 
cases, respectively. 

ties intermediate between those of the Si-, 0-, and OO- 
terminated cases, we enlarge the supercell by a factor 
of two along the y-axis (along the horizontal direction 
of Figs. 1, 2, and 6). For the lx2-supercells, three new 
surfaces, in addition to those already considered above, 
are now possible. These are: a (Si,0)-terminated sur- 
face which consists of neighboring Si- and O-terminated 
surfaces; a (Si,00)-terminated surface which consists 
of neighboring Si- and OO-terminated surfaces; and a 
(0,00)-terminated surface which consists of neighboring 
O- and OO-terminated surfaces. [Note that the (Si,00)- 
terminated case possesses the same oxygen density as the 
O-terminated case.] We find that the (Si,00)-terminated 
substrate has energy which is higher than that of the O- 
terminated surface by ~ 1.4 eV per Si atom at the surface 
(per 25 A^). Therefore, we restrict our focus to investi- 
gation of the (Si,0)- and (0,00)-terminated cases, cal- 
culating their relaxed structures and adhesion energies. 

Figure 7 depicts the adhesion energy as a function of 
the oxygen density at the interface. Remarkably, the ad- 
hesion energy is not a monotonically-increasing function 
of oxygen density; at low oxygen density, increasing the 
number of O atoms a small amount does not necessarily 
lead to stronger adhesion. However, beyond the oxygen 
density of 4 nm~^ where the 00-termination will begin 
to appear appreciably, the adhesion energy increases, and 
does so almost linearly. This linear dependence implies 
that each 00 surface unit contributes to the adhesion 
fairly independently, and thus the locally OO-terminated 
regions may not necessarily have to form large domains 
in order to work as an effective "glue" , although it is clear 
that a critical density of such locally OO-terminated re- 
gions is required for strong adhesion. 
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FIG. 8. Relaxed structures of the OHOH-terminated sur- 
face and interface. 

3. Hydroxylated surfaces 

We now address a frequently-cited reason for the weak 
adhesion at Cu/a-Si02 interfaces, namely the existence 
of hydroxyl groups on the Si02 surface.^ R is well known 
that hydroxyls are likely to form on the Si02 surface un- 
der "wet" deposition conditions; to examine the impact 
of hydroxyls on adhesion from first-principles, we add two 
hydrogens to the 1x1 OO-terminated surface, and relax 
all of the atoms. The result is that the cohesive energy 
of OHOH-terminated substrate, in which a Si atom is 
terminated by two hydroxyls, is lower than that of OO- 
terminated substrate by ~ 11.2 eV per surface Si atom. 
This energy difference is much larger than the dissoci- 
ation energy of H2 molecule (~ 4.5 eV per molecule), 
suggesting the high stability of the OHOH-terminated 
substrate. If we now place the Cu layers on top of the 
hydroxylated surface and relax the interface (Fig. 8), we 
find that the deposition of Cu layers does not significantly 
affect the structure at the OHOH-terminated surface, im- 
plying little interaction; the hydrogen present at the sur- 
face leaves the surface neutral and inert. This tendency 
is clearly reflected in the calculated adhesion energy, and 
according to our calculations it is just 0.331 J/m^, less 
than one tenth of the adhesion energy observed in the 
OO-terminated case. 

However, as shown in Table II, the adhesion ener- 
gies for all other cases studied here (i.e. the 00-, O- 
, and Si-terminated cases) seem to be of a magnitude 
which should produce adhesion, especially the oxygen- 
rich OO-terminated case. It is also noteworthy that the 
OO-terminated substrate, which has the largest adhesion 
energy among the three, looks actually the most stable 
substrate. These facts imply that should dehydroxylation 
of the substrate surfaces be successfully achieved (e.g. by 
annealing or particle bombardment of a growing film sur- 
face), robust adhesion is entirely possible, at least from 
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the point of view of chemical bonding. 



V. CONCLUSIONS 

We have performed a first-principles study of the ad- 
hesive properties of atomically-sharp Cu/Si02 interfaces. 
As a model a-Si02 substrate, we used a crystalline poly- 
morph of Si02, a— cristobalite, and investigated its Si- 
, 0-, and 00-terminated (001) surfaces in detail. For 
interfaces between Cu and 00-terminated surfaces, a 
substantial rearrangement of oxygen positions relative 
to the free surface is observed, suggesting the formation 
of strong Cu-0 bonds. Analysis of the local density of 
states at the interface showed that Cu-0 interfacial bonds 
are composed of Cu—d and 0—p states. The computed 
adhesion energy also exhibited a tendency toward much 
stronger adhesion in the 00-terminated interface than 
in the O- and Si-terminated interfaces: Substrate sur- 
faces with high oxygen content were found suitable for 
adhesion. 

The adhesion energy is found to be very insensi- 
tive to the position of the Cu layers if the termination 
type is unchanged, and this observation suggests that 
a— cristobalite may serve as a good starting model of 
a-Si02 substrates. The detailed dependence of the adhe- 
sion energy on oxygen density at the substrate surface has 
been investigated, and it shows that the surfaces where 
the O- and 00-terminations coexist can also lead to rel- 
atively strong adhesion. 

The possible existence of hydroxyl groups at the sub- 
strate surface is thought to be the main cause for the 
weak adhesion observed in experiment. However, if these 
hydroxyls are removed beforehand and the oxygen den- 
sity at the surface is increased (for example, by upping 
the oxygen partial pressure), Cu layers are predicted 
to adhere reasonably well to Si02. In addition to the 
possibility of direct comparison to future measurements 
on crystalline Cu/Q:-cristobalite(001) interfaces, these re- 
sults should serve as a starting point for which to study 
more complicated interfacial geometries, benchmarks for 
future investigations, and a baseline for experiments at- 
tempting to elucidate complicated phenomena affecting 
adhesion at Cu/glass interfaces. 
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